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Abstract 

The Heston model is a popular stock price model with stochastic 
volatility that has found numerous applications in practice. In the 
present paper, we study the Riemannian distance function associ- 
ated with the Heston model and obtain explicit formulas for this 
function using geometrical and analytical methods. Geometrical ap- 
proach is based on the study of the Heston geodesies, while the 
analytical approach exploits the links between the Heston distance 
function and the sub-Riemannian distance function in the Grushin 
plane. For the Grushin plane, we establish an explicit formula for 
the Legendre-Fenchel transform of the limiting cumulant generat- 
ing function and prove a partial large deviation principle that is true 
only inside a special set. 

1. Introduction 

There are two main protagonists in this paper: the Riemannian 
manifold associated with the Heston model of the stock price, and 
the Grushin plane, which is one of the best-known examples of a 
sub-Riemannian manifold. The present paper focuses on the Heston 
Riemannian distance and the Grushin sub-Riemannian distance and 
provides explicit formulas for these. The Heston distance and the 
Grushin distance are intimately related, and various facts concerning 
these distances can be easily transplanted from one setting into the 
other. 

We will next briefly describe the main results obtained in this pa- 
per. Theorems 1 and 3 below contain explicit formulas for the Heston 
distance. The formulas in Theorem 1 are established using geomet- 
rical methods, while the proof of the distance formula in Theorem 3 
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uses certain links between the Heston and the Grushin distances and 
is more analytical. In the proof of Theorem 3, we compute and study 
the limiting cumulant generating function A for the Grushin plane 
and the Legendre-Fenchel transform A* of the function A. One of the 
main results in the present paper is a partial large deviation principle 
for the Grushin plane (see Theorem 22). The word "partial" is used 
in the previous sentence because in the case of the Grushin plane 
the large deviation principle with A* as a rate function holds only 
inside a special subset of x R^. We would also like to bring the 
reader's attention to the results concerning certain qualitative prop- 
erties of the transcendental equations whose solution is involved in 
determining the Heston Riemannian distance function. These con- 
vexity and monotonicity properties established in Lemmas 8, 9, and 
10 ensure that the equations can be efficiently and rapidly solved 
by Newton's method or a bisection method. We also show in the 
present paper that it is crucial to distinguish two different regimes 
(the near and the far point regime) in the geometrical and analytical 
approaches to the Heston distance, each regime requiring it's own 
analysis (see Theorems 1 and 3 and their proofs). 

Let us expand on the financial motivations for considering the 
Heston Riemannian distance function. The Heston model is one 
of the most popular stock price models with stochastic volatility. 
This model was introduced in [ ]. More information on the Hes- 
ton model and stochastic volatility models can be found in [15, 16, 
20, 31]. The stock price process S and the variance process V in the 
Heston model satisfy the following system of stochastic differential 
equations: 

f dSt = }iStdt + y/YtStdWt , . 

\ dVt = {a- bVt)dt + c^/VtdZt, ^ ' 

where fl > 0, > 0, c > 0. In (1), W and Z are correlated standard 
Brownian motions such that d{W, Z)t = pdt with p G ( — 1, 1). 

Recently closely related models, the local-volatility Heston mod- 
els, given by dSt = jiStdt + y/VtO'{St)StdWt for an appropriate func- 
tion a, have become objects of widespread interest among practi- 
tioners. Practitioners seek to come up with accurate approximations 
to the Black-Scholes implied volatility in such models and there is a 
considerable literature in this direction. In one of the approaches to 
this problem, initiated for another class of stochastic volatility mod- 
els in [17], and in the Heston case by [ ], [13], [14] a key element 
in determining the implied volatility is the Riemannian distance to 
a line S = K in the SV-plane. We plan to address applications of 
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the results obtained in the present paper to the local volatility Hes- 
ton models in future publications. This will include consideration of 
heat kernel expansions, implied volatility expansions, and pricing of 
exotic options in such models. 

Let us consider the following uncorrelated Heston model: 

dSt = StVVtdWt 

dVt = {a- bVt)dt + VVtdZt, ^ ' 

where a > 0, b > 0, and W and Z are independent standard Brown- 
ian motions. Denote by X the log-price process defined by X = log S. 
Then the model in (2) transforms as follows: 

dXt = -\Vtdt + ^tdWt 



dVt = {a- bVt)dt + ^idZt. 
The state space for the process (X, V) is the closed half-plane 
U = {(x,z;) G ]R2 . ^ > q|_ 

We will denote the initial condition for the process (X, V) by {xq, vq). 

The Riemannian metric form associated with the Heston model is 
defined on the interior Ti° oi the closed half -plane Ti as follows: 

ds^ = (^dx^ + dv^^ . (4) 

The open half-plane Ti° equipped with the metric defined in (4) is 
called the Heston manifold. The form in (4) generates the Riemann- 
ian distance dn onTi. More examples of Riemannian distances aris- 
ing in finance can be found in [ ] (see also [ ]). 

In this paper, we discuss various explicit formulas for the Heston 
distance dn- It is worth mentioning that the following two-sided 
estimate for dn is known (see [21], Proposition 4.3.2): 

D{xo,vq,xi,vi) < dH{{xo,vo),{xi,vi)) < 12D(xo,^0/^i/^i) (5) 
for all {xq,Vq) G T-L and {xi,Vi) G H, where 



We will next briefly explain how to obtain explicit formulas for 
the distance function dn in the general correlated Heston model de- 
scribed in (1) from similar formulas for the distance function dn in 
the uncorrelated Heston model with the vol of vol coefficient equal 
to one considered in the present paper. It is known that the principal 
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part of the generator of the diffusion (X, V) with X = log S in model 
(1) is given by 

2 c-^ c2— 
dx^ dxdv 9z;2 

Let us show how to reduce this to the case of an uncorrelated Heston 
model, in which the metric is in the standard form presented in (4). 
First make the change of time t = c^t which reduces the principal 
part of the generator of the diffusion to 

It is not hard to see that under such a change of time, the distance 
function rf^ is multiplied by the constant c. Next, we use the diffeo- 
morphism 

c p 
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V. 



After this change of variables the principal part of the new diffusion 
operator is in the standard form 

+ 




9x2 3^2 

Taking into account the reasoning above, one can prove that 

dHiixo,vo),(xi,vi)) 

pVQ \ ( cxi- pvi \ \ 

Formula (7) shows how to adapt the distance formulas obtained in 
this paper to the case where the Heston model is correlated and 
given by (1). 

Let us note that the drift terms in (3) do not affect the Heston dis- 
tance. On the other hand, the Heston transition density {x,v) ^ 
p|^((xO/^o)/ (^/^)) associated with the process (X, is influenced 
by the drift terms. Using the definition of the Riemannian distance 
and (4), we see that the distance satisfies the following conditions: 

dniixo, vq),{xi,vi)) = dH{{Q,vo),{xi - xq,vi)) (8) 

and 

dH{{oiXo,ciVQ),{ciXi,ciVi)) = y/ocdH{{xo,vo), {Xi,Vi)) (9) 
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for all a; > 0. 

It will be assumed throughout the paper that Xq 7^ Xi. In the case 
where xq = x\, the geodesies joining the points {xq,vq) and {x\,V\) 
is a vertical line with length 

/■max(Do,Ui) 1 j j — \ 



It follows from the previous formula that the x-axis is at finite dis- 
tance from any point in "H", and hence the x-axis, being part of the 
boundary of the Heston manifold, is "at infinity". Therefore, the He- 
ston manifold is not complete, and we can not apply the Hopf-Rinow 
theorem to establish the existence of a length-minimizing geodesic 
joining two points in H. Note that one difficulty in establishing such 
a result is proving that the length minimizing curve joining any two 
points is a true geodesic and not the union of broken geodesies (on 
the other hand for the metric ds^ = v{dx^ + dy^), Robert Bryant has 
communicated to us that only a subset of points in the upper half 
space can be joined by a non-broken geodesic). 

It is interesting that for the metric defined by (4) the existence and 
uniqueness result for the length-minimizing geodesies has been es- 
sentially known for at least one century. We will next provide more 
information. 

Let us consider the following metric form: 



that is intimately related to the Heston metric. It is known that the 
length-minimizing geodesies for this metric are dilations and shifts 
of the standard cycloid given by 



This was established by O. Bolza in 1904 (see [2], see also [18], Propo- 
sition 1.2.1). Bolza proved that for any two points C = {xq,vq) and 
D = {xi,Vi) in % with xq 7^ X\ there exists exacty one dilated and 
shifted standard cycloid joining C and D and such that no cusp of 
the cycloid lies between C and D. Bolza states that the previous re- 
sult was included without proof in unpublished lectures of K. T. W. 




and hence 



^^h((^0/^'o)/(^0/^'i)) = 2(Wmax(z;o,Ui) 






(10) 
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Weierstrass (1882). A special case was handled by Johann and Jacob 
Bernoulli at the very end of the 17th century. 

It follows from the above-mentioned result of Bolza that the Hes- 
ton geodesies can be obtained from the cycloid curve 

1"=^'"'^"* (11) 
1 u = 1 — cos 

by translation x x + b and dilation x cx, v cv, or are vertical 
lines. It is also clear that Bolza's existence and uniqueness theorem 
holds for the Heston manifold. 



2. The Heston distance 

For every fixed C G R such that C 7^ 0, define a function on the 
interval [0, C~^] by 



Now fix vq>Q and Ui > and consider a function of the variable C 
given by 

-F{vQ,vr,C) =F{vi,C) -F{vq,C). (13) 
It is assumed in (13) that the variable C satisfies the condition 

1 
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< |C| < min yv^ ^ ,v-^ 

Since lim Tivn, Vi, C) = (use I'Hopital's rule twice), we can extend 

c^o 

the function F continuously to an odd function on the interval 

_i _i _i _i 

1 = [-min(z;o ^,v^ ^),rmn{vQ ^,v^ 2)] (14) 

by putting F(uo, Ui,0) = 0. 

We will next prove a theorem that provides explicit formulas for 
the Heston distance. Note that there are two expressions for the He- 
ston distance in Theorem 1, depending on the location of the points 
{xq,vq) and {xi,vi) in the Heston half-plane. 

Theorem 1. (i) Suppose the points {xq,vq) and {x\,V\) in the Heston 
half-plane satisfy the following condition: 

1 

1^:1 -xol <F{min{vQ,vi),max{vQ,vi),—= ), (15) 

^Jm.ax[VQ,Vl) 



where F is given by (13). Then 
dH{{xo,vo),{xi,vi)) 

arcsin(C* ^/max{vo,Vi)) — arcsin(C* ^ymin(vo,Vl)) 



(16) 



C* 

where C* = C*{{xo,Vo), {xi,vi)) is the unique solution to the transcen 
dental equation 



a/uov1 — C^vq axcsmiC^/vQ] 
^1 - ^0 = — 



C C2 



Vw/l — C^Vi arcsin(Ci/ui' 



C C2 
on the interval I defined by (14). 



(17) 



(ii) Suppose the points {xq, vq) and (xi, z^i) in the Heston half-plane satisfy 
the following condition: 

1 

\xi-Xo\ > F{min{vo,Vi),max{vo,Vi), ). (18) 

y/max{vo,Vi) 

Then 

^^((^^o, ^^o), [Xi, vi)) = 2 ^^-^^ ^ i^, (19) 

where C* = C*{{xo,vq), {xi,vi)) is the unique solution to the transcen- 
dental equation 



vq^/I — C^vq arccos(Cy^ 



xi - XqI = h 



C C2 



V^yl -C2z;i arccos(Cv^z;i) 
^ C C2 ^ ^ 
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on the interval [0,min(z;Q ^,v-y ^)]. 

Proof. We will first prove part (i). It is not hard to see that with no 
loss of generality we may assume that xq < x\ and vq < vi. Let c 
be the dilation coefficient corresponding to the points A = {xo,Vo) 
and B = {xi,vi) in Bolza's description of the geodesies. Then we 
have c > 0. Define C > from the equality c = 2C^, and consider 
the geodesic, connecting A with B, under the scaled parametrization 
s Cs. 

Let us first assume that the point A and B are both to the left of the 
apex of the arc of the geodesic passing through them. Then it is not 
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hard to see that the components s ^ x{s) and s i— > v{s) of the arc of 
the geodesic through A and B satisfy 

X = Cv 

V = v^Vl - C^v, ^^^^ 

where the derivative in the system above is taken with respect to the 
arclength parameter and where v G [0, C~^] . It is also clear that x is 
a function of v for v G [0, C~^] and we have 

dx ^ 

dv ^l-C^v 

Note that v = is the second component of the apex. Therefore 

[vQ,vi] C [0,C-2]. 
It is easy to see that 

f- = ■ (23) 

It follows from (22) and (23) that f (u, C) = x{v) -\- a, where a is some 
constant. Plugging v = Q into the previous equality, we get a = 
b, where b is the shift parameter in the description of the geodesic 
passing through A and B. Hence, 

F{v,C) = x{v) + b. (24) 
It follows from (24) that C satisfies the following condition: 

xi-xo = F(vi, C) - f (570, C) = F{vo, vi, C). 

Moreover, the assumptions formulated above can be formulated as 
follows: 

xi-xo <F ( Vo, vi, 



It is not hard to see using (22) that 



Wo 




^h(a,b)= Wb:- +1 



^ dv 1 dv 



\fv Jvo yj\ - C^V Vv 



^^Lrcsm{C^yvl) — arccos(C^z;o) 

c ■ 

This establishes formula (16), and completes the proof of part (i) of 
Theorem 1 in the case where the points A and B are both to the left of 
the apex of the arc of the geodesic passing through them. The proof 
of part (i) in the case where A and B are to the right of the apex is 
similar. 
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Next suppose that one of the points is to the left of the apex, while 
the other one is to its right. This case is a combination of the previous 
two. With no loss of generality, we may assume that (xq, vq) is to the 
left of the apex and {xi,vi) is to the right of the apex. This happens 
if and only if condition (18) holds. It is not hard to see that under the 
restriction imposed above, we need to sum two contributions, one 
going from {xq,Vq) to the apex and the other going from the apex to 
(xi, Vi). Since the u-component of the apex equals C~^, we obtain 

_ y/vo^/l - C^vq 

^apex •^0 — Q 

_ y/VQ\/l - C^VQ 

c 

To this we must add 

Xl Xapex — ^ 

c 

so, in total we obtain condition (20) for C. In the same way we must 
add the corresponding distance formulas to get formula (19). 
This completes the proof of Theorem 1. 

We will say that the points (xo,i/o) ^rid are C-close pro- 

vided that the inequality in (15) holds. Similarly, if the inequality in 
(18) holds, then we will say that the points are C-far. The equations 
in (17) and (20) will be called the C-equations, while the formulas in 

(16) and (19) will be called the C-formulas. 

Remark 2. Two different formulas appear in Theorem 1 because cer- 
tain subtleties which underlie the geometry of the cycloid have to be 
dealt with. Note that formula (16) was suggested as the Heston dis- 
tance formula in the book [ ] by P. Henry-Labordere (see formula 
(6.66) in [ ]). However, formula (16) holds only in the close-point 
regime (see part (i) of Theorem 1) and has to be replaced by formula 

(17) in the far-point regime. The presence of two different regimes 
was not taken into account in [20]. 

It follows from part (ii) of Theorem 1 that 

dH{{xQ,0), (xi,0)) = 2^ n\xi - xo\. (25) 

The previous formula describes the Heston distance between any 
two points on the boundary of the Heston half-plane. 



arcsin(Ci/yo) n 
C2 ^2C2 
arccos(Cy^z^) 

arcsin(Ci/yi) n 
C2 ^2C2 
arccos(C1y^^i") 
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3. The toy Heston model and the Grushin model 

Consider the following stochastic model: 

dGt = HtdWt 
dHf = ^dZt- 



(26) 



We will call the model described by (26) the Grushin model because 
the Laplace operator associated with it is a special Grushin operator 
given by 

L = l|i + y2|i. (27) 
4 91/2 -^9x2 

The heat kernel for the operator L will be denoted by ( (xq, i/o), (xi, i/i ) ) 
where ^ > 0, {xq, j/q) £ IR^/ and (xi, i/i) G IR^, and the sub-Riemannian 
distance on IR^ Carnot-Caratheodory disrance), corresponding to 
the Grushin model will be denoted by rfg. The plane IR^ equipped 
with the distance dg is called the Grushin plane. The heat kernel 
for the Grushin plane satisfies the following partial differential 
equation : 

-df - + 2^ ^^^^ 

with the initial condition given by 

((^o,yo),(:^i,yi)) = ^xo{xi)5y^{yi) 

More information on the geometry of the Grushin plane can be found 
in [6, 7, 8, 9, 27, 28]. Stochastic methods which are used in the study 
of Grushin type structures are discussed in | |. 

Suppose the process (G, H) is the solution to the system in (26) 
with the initial conditions gQ and Hq, respectively. Then the process 
(X, y), where X = G and Y = H^, solves the following system of 
stochastic differential equations: 

dXt = ^/YtdWt 



dYt = \dt + ^/YtdZt ^^^^ 

with initial conditions go and h^. In (29), the processes W and Z are 
new standard Brownian motions defined by Wf = sign (Hf)rfWt and 
Zt = siga(Ht)dZt. We will call the stochastic model described by 
(29) the toy Heston model. In this section, the Heston distance will 
be analyzed by using the following formula relating the Heston and 
the Grushin distances: 

dH{{xo,Vo),{xi,Vi)) = dG{{X0,^/v^),{Xl,^/vi)), (30) 
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for all points A = {xq,Vq) G % and B = {xi,Vi) G %. It is not hard 
to prove equality (30) when A and B belong to and then extend 
the equality to % by continuity. The proof for 1-1° is based on the 
fact that the length minimizing Heston and Grushin geodesies for 
A G H° and B ^ Ti° are entirely contained in We leave filling in 
the details as an exercise for the reader. It follows from (30) and (5) 
that 

D{xo,vo,Xi,Vi) < dciixo, y/v^), {xi, y/vl)) < 12D(xo,yo/^i/^'i) 

(31) 

for all {xo,vo) G H and {xi,vi) G H, where the function D is given 
by (6). 

We will next formulate a statement which provides an alternative 
formula for the Heston distance. 

Theorem 3. For any two points {xq,vo) G H and {xi,vi) G H such that 
at least one of them is not on the boundary, the following formula holds: 



duiixo, vo),{xi,vi)) 



sin 



vi + vo- 2^viVo cos I - I , (32) 



where S = S{{xq, vq), (xi, Vi)) is the unique solution to the equation 



[vi + Vq) {5 - sin(^)) - [d cos (^|) - 2 sin 



2sin2 ( I 



Xi - Xq, 



(33) 



satisfying the condition —In < S < 2n. 



Remark 4. In [27, 28], M. Paulat established a formula for the sub- 
Riemannian distance in a slightly different Grushin model given by: 



dGt 
dHt 



HtdWt 
dZt. 



Paulat's formula is equivalent to (32) (one formula can be obtained 
from the other using (30)). The ideas used in the proof of formula (32) 
are completely different from those employed in [27, 28]. Paulat an- 
alyzes sub-Riemmanian geodesies in his proof, while the techniques 
used in the present paper are more analytical. In addition, the proof 
of Theorem 3 contains several new results, e.g., a partial large de- 
vation principle for the Grushin model. We would like to thank M. 
Paulat for sending us his dissertation [28]. 
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Note that there are two distance formulas in Theorem 1 (in the 
close point regime and in the far point regime), while Theorem 3 
contains only one distance formula. An interesting fact is that in the 
(5-environment, there is a special two-set partition oiTi x Ti hidden 
in the background. 

Definition 5. We will say that the points {xq,vq) and {xi,vi) are 5-close 
■provided that 

1^1 - ^o| < 2" (^1 + ^'o) + 27^1^0/ (34) 

and 5-far if 

1^1 - ^o| > y (^1 + ^o) + 2V^iz;o. (35) 

Remark 6. In terms of the parameter d, the description of the close- 
point (^-regime and the far-point ^-regime is < |^| < n and n < 
\S\ < 271, respectively. 

Remark 7. It follows from Theorem 1 or Theorem 3 that 

du {{xo,vo), (xi,ui)) = dn {(xi^vo), (xo,ui)) 

and 

du {{xo,vo), {xi,vi)) = dn ((xo,ui), {xi,vo)) . 

The reasons, why the two regimes in Definition 5 are introduced, 
are rather subtle. It will be shown in the next sections that the close- 
point ^-regime describes those pairs of points, for which formula (32) 
can be obtained by analytical methods. The far point ^-regime is a 
proper part of the far point C-regime, and formula (32) in the far 
point regime can be established using formula (19) (see Lemma 14 
below). 

We will next derive formula (25) from Theorem 3. With no loss of 
generality we can assume xq < x\. Let e > 0, and take = 0, Ui = £. 
Then Theorem 3 implies that 

-1 



— ^ 

5 



where 



Therefore, 



dH((^0/0),(xi,£)) = 5 



£{S — sinS) 
2sin2| 



dH{{xo,0),{xi,e)) 



sm I ^ 



Xl - Xq. 



2{xi-J 



Sir 
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It is not hard to see that S ^ 2n as e 0. Taking the limit as e 4- in 
formula (36), we obtain formula (25). 

4. Solvability and Convexity 

In this section we discuss the unique solvability of the C-equations 
and of the cJ-equation. Let us start with the C-close point regime. It 
is clear from the definition of the function F in (13) that to study the 
unique solvability of equation (17), it suffices to assume Xi > xq and 
Vi > vq. Then the equation becomes 

F(vo,Vi,C) = xi- Xq, (37) 

and we have to solve it on the interval < C < z;^ ^. Note that in the 

C-close point regime we have Xi — xq < F{vo, V\, ^), and the func- 

tion f maps the interval [0,u-^ ^) onto the interval [0,f (uq/^I/^i ^)) 
(the monotonicity of the function F follows fom the next lemma). 

Lemma 8. For fixed V\ > vq, the function C i— F{vo,Vi,C) is strictly 

increasing on the interval [0, ^ ) and convex on the interval (0, z^^ ^ ) . 

Proof. Consider the function 



arcsm{Cy/v) y/vVl - C^v 
h{v,C) = ^2 ^ . 

By definition, the function F satisfies 

F{vo,vi,C) = Fi{viX)-Fi{voX). (38) 



Therefore 



d^F(vo,v,,C) d^Fi a2fi 

— ac^ = ac2(^i'C)-^(.o,c) 

By differentiating the function Fi, we get 

dFi _ 2 arcsin(C-\/u) ^ y/v 



and 



dC C3 CVl - C^v 



3% ^ 2V^(-3 + 4Cy) ^ 6arcsin(Cv^) ^^^^ 

C^(1-CW ^ 
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We also have 



^'f.;,c)^^!E2i^. ,42) 



This is quite remarkable because if we stop at the second derivative 
we have the more complicated expression given in (41). Now, it is 
not hard to see, using (39) and (42), that 



ac2 



> (43) 



for all C G ^0, v-^^j, and hence the convexity statement in Lemma 
8 holds. 

We will next prove that the function F is increasing. Using (38) 
and (40), and making tedious but straightforward computations, we 
obtain 

|=:lim|(.o..,.C) = ?{.f-4)>0. (44) 

Now the fact that the function F is increasing follows from (43), (44), 
and the equality 

— (z;o,ui,C) = ^[vo,vi,u)du + —{vo,vi,0). 
This completes the proof of Lemma 8. 

It follows from Lemma 8 that equation (37) is uniquely solvable 
for all pairs of points in the Heston half-plane which are C-close. 

Our next goal is to prove a similar result for any pair of points 
{xo,Vo) and (xi,Vi) in the Heston half-plane which are C-far. With 
no loss of generality we assume xq < ^i arid z^o < Vi. Recall that the 
C-far point regime is described by the following inequality: 

xi - xo> F (vo,vi, ^ 



The C-equation in the far point regime is as follows: 

F(uo, vx, C) = xi - xo, (45) 

where 



~ v^^/l - C^VQ arccos(Cyz;o) 
F{vq,Vx,C) = ^ h 



C C2 
^ ^y/l - C^vx ^ arccos(CVz^) ^^^^ 
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_ 1 

and we are looking for the solution C* satisfying < C* < v.^^. 

1 ' 



The function f is decreasing on the interval |^0, v.^^ j (see Lemma 
9 below), and maps it onto the infinite interval (^Vq, Vi, j / 



Indeed, it follows from the definition of F that \ivs\F{vr),v-\,C) = oo. 

c;o 

We also have 



lim F(vo, V\, C) = \/vnV\^ /l — — + z^i arccos 




= V^^Wl h -^1 -^larcsmW— = f uo,z;i,— ^ 

V ^1 2 V ^1 V \/W 

It is clear from the previous discussion that equation (45) is unquely 
solvable in the C-far point regime and the solution C* satisfies < 

C* < v~^. 

Lemma 9. Let vq <v\. Then the function F defined by (46) is decreasing 

on the interval ^0, ^ ^ • locally convex near the point C = and 

_i 

locally concave near the point C = v-^^. 
Proof. Put 



F, („, c) = N^^V^ + (47) 

Then we have 

F{vo,ViX) = Fi{voX) + F{vi,C). (48) 
By differentiating the function F twice and simplifying, we obtain 



2^/v ^ 2 arccos(Cv^) 



and 

B^fi 6^/v — 8v2C^ ^ 6 arccos (Ci/z;) 

C3(l-C2z;)i 



(49) 



It is clear from (49) that the function C i— > fx(z;, C) decreases. 
Moreover, formula (50) shows that 
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and 

1 dC^^^'^^ = -oo. 

Analyzing the previous equalities and taking into account (48), we 
see that Lemma 9 holds. 

Let us next consider the equation in (33). We fix xq G R, xi G R, 
^0 ^ 0/ ^1 ^ ^rid assume that at least one of the numbers z^o and vi 
is different from zero. Denote the function on the left-hand side of 
(33) by /(^). Then the equation in (33) can be rewritten as follows: 

f{S) = Xi- Xq. (51) 



We have 



where 



and 



f{S)=f{S,Vo,v^) = ^A{S)B{S) (52) 
A{S) = A{S,Vo,v-,) = sin-2 (^0 (53) 



B{S) = B{S,Vo,vi) = {vi + vo) {S-sm{S)) 

-2,/^^scos^^-^ -ismf^^y (54) 



The value of the function / at ^ = is given by 

fiO) = UimA{3)B{S)=0. (55) 

The function / is continuous and odd on the interval (— 2/1,2 tt). 
In order to prove that the equation in (33) is uniquely solvable on 
{—In, In), it suffices to assume that Xi > xq and look for the unique 
solution belonging to the interval (0, 2n). We will next show that the 
function / is positive on the interval {0,2n). Indeed, the functions 

3 5 — sin(^) and 3 3 cos — 2 sin are equal to zero at 

3 = 0. Moreover, the former function is positive on {0,2n), while 
the latter one is decreasing (differentiate!) and hence negative on 
{0,2n). It follows from (54) that B{S) > for all S G {0,2n). This 
shows that the function / is positive on the interval (0, 2n). 
We also have f{n) = j[n {vi + Vq) + A^viVq]. 

Lemma 10. The function f is strictly increasing and convex on the interval 
(0,27r). Moreover, it maps (0,27r) onto (0, oo). 
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Remark 11. Lemma 10 implies that there exists the unique solution 
to the equation in (33) belonging to the interval (—271, 2k). 

Proof. For all 5 G (0, 27r), we have 

A'(^) = -sin-3(^0cos(0 (56) 

and 

2sin4 U\ 

B'{5) = (ui + vq) (1 - cos(^)) + v^'Jsin (^0 , (58) 

and 

B"{5) = {vi + vo) sm{S) + (sin (^0 + cos (0 ^ . (59) 

Next, using the product rule, THopital's rule, and the formulas above, 
we see that 

lim/((5) = Uvi+ Vo) + ^ (60) 

Now, (55), (60), and the fact that the function / is odd on (—In, In) 
imply that / is differentiable at S = 0. 
Our goal is to show that 

f"{^)>0 for ^G(0,27r). (61) 

It is not hard to see that if (61) holds, then the function / is increasing 
on (0,27r). Indeed (61) implies that for all < e < (5 < 2n, 

f'iS)=f'ie) + lj"iu)du>f'{e). 

Therefore, (60) shows that the derivative of the function / is posi- 
tive on (0,27r). Now the continuity of / on [0,27r) implies that the 
function / is increasing on [0, 2n). 

Our next goal is to prove (61). We have 

(AB)" = A"B + 2A'B' + AB" . (62) 

Using (53), (54), (56), (57), (58), and (59), we obtain 

(2 + cos((5)) 



A"B 



2sinM| 



[vi + Vq) {5 - sm.{5)) -2./z^[t^ l^^cos |^0 -2 sin 
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2A'B' 
and 

AB 



2 cos 






{vi + Vq) (1 - cos(^)) + y'vivod sin 






sin^ 







(ui + ^o) sin(^) + y'z^Iz^ (^sin + jcJcos (^1 

.2 f s 



sin I 2 



It is not hard to see using (62) that (61) holds if the function 

/(cJ) = (2 + cos(cJ)) 



4 sin 



^1 + ^o) ~ sin(^)) — 2y/viVQ ^^cos — 2 sin 

cos (t^i +uo) (1 - cos(^)) + V^i^o^ sin 

+ 2 sin^ (yi + z;o) sin(^) + ^/v^ (^sm (^0 + cos (^0 

is positive. The previous function can be represented as follows: 



where 
and 



Ci ((5) = 2(5 + ^ cos(^) - 3 sin((5) 



C2(.5) = 6sin(0 ( 



1 + cos^ ( - 



S COS - 5 + cos - 



(63) 



Now it is clear that it suffices to show that Ci{S) > and C2{S) > 
for all (5 G {0,2n). 
We have 

Ci(o) = c;(o) = ci'(o) = o 

and 

C'"(^) = (5sin((5) > 

for all ^ G (0,27r). Therefore Ci(^) > on the interval (0,27r). 

In order to prove that C2{S) > for all < S < 2n, we rea- 
son as follows. First note that (63) implies the inequality C2{S) > 
for all 3 G [n,2n). Therefore, we can restrict ourselves to the case 
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where < ^ < tt. The following estimates, which can be easily de- 
rived using Taylor expansions of the sine and cosine functions, will 
be needed in the proof below: 

> - 



sm 



48' 



^«^'-'<^-^ + 384' 



and 



2^ 

sin 5 < 5 



for all (5 G (0,7r). 

Using (64) and (65), we obtain 



> 5 
= 5 

> 5 



6 sin 

'4 



1 + cos^ 



2 
2 



+ 2cos^ 
64 



1 + cos^ 

5_ 
2 



5^ 5^ 
6"^ 120 



5 cos 
, 5^ 



5 + cos 



8 ^ 384 



5 + cos^ 



384 



5 + cos 



2 sin 



It will be shown next that the function 

h{5) = — - — -2sin^ - 
^ ^ 2 64 V2 



is positive on the interval (0, zr). Indeed, h{Q) =0 and 
Now using (66) we get 



5^ 

5 sin 0. 

16 



120 



> 



(64) 
(65) 

(66) 



(67) 



since 5 G (0, n). It follows that C2{6) > for < ^ < tt. Finally, it 
is easy to check, using the definition of the function /, that / maps 
(0,27r) onto (0,oo). 
This completes the proof of Lemma 10. 

Lemma 12. Suppose the points {xo,Vo) G Ti and {xi,vi) G are 5-far. 
Then they are C-far. 
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Proof. It suffices to prove the lemma in the case where Vi > 0, 
^^0 > 0, and Vi > vq. The assumption in the formulation of Lemma 
12 means that inequality (35) holds. To prove that the points are C- 
far, we have to establish estimate (18). For vi > vq, this estimate is 
as follows: 



Ixi — xqI > + J viVq — — Ui arcsin . /— . (68) 
2 V " y Vi 



It is easy to see that condition (35) implies condition (68). 
This concludes the proof of Lemma 12. 
Let us define the following function: 

sin — 

It is assumed in (69) that u G {7i,l7i), Vq > 0, and Vi > 0. We also 
exclude the case where vq = vi = 0. It is not hard to see that for 
fixed Vq and Vi, the function u i— > f{u, vo,vi) is strictly decreasing, 
continuous, and 



/: (71,271)^ (o, J—- ) 



(70) 



(the mapping is onto). 

Suppose the points (xq/^o) ^ ^ and {xi,vi) G H are S-far. By 
Lemma 12, these points are also C-far. It is natural to ask wheter 
there exists a relation between the numbers S* and C*, correspond- 
ing to the given points. The next statement answers the previous 
question. 

Lemma 13. If {xq,vo) G 71 and {x\,Vi) G % are 5-far, then 

C*{{xq,vq),{xi,vi)) = / (|(5((xo/^'o)/(^i/^^i))|/^'0/^'i) / (71) 

where the function f is defined by (69). 

Proof. Since the points are ^-far, we have n < \S\ < In. With no 
loss of generality, we can assume that n < 5 < In, xq < ^1/ arid 
< Note that condition (70) implies that all the expressions 
appearing in the proof of Lemma 13 are real numbers. 
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— , and recall that the number 5 is the unique 



Let us denote a = 
solution to equation (33). Then 

Xi - Xo 

(1 + (x^) {S - sin((5)) - 2(x (^Scos (^^S^ - 2 sin (^^S 



2sin2 ( 1(5 



(72) 



Put 



C = v^/ {5{{xo,vo),{xi,vi)),vo,vi) 
sin ^2^^ 



l + a;2-2acos (^^^^ 
It is not hard to see that equality (71) holds if and only if 

Xl — Xq 



(73) 



ay 1 — a^C^ + v 1 — ^ arccos C + arccos(aC) 



C 



C2 



(74) 



where C is defined by (73). Using the addition formula for the in- 
verse cosines, we see that equality (74) is equivalent to the following 
equality: 



Xl-Xo I 1 2 , / 1 1 

C2 VC2 



^^1 



+ 



arccos (a:C2 — — C^\/l — a'^C'^) 



C2 



Since ^ solves equation (33), we have 

Xi - Xq 



[1 + a.^){S - sin (5) + 2oi (2 sin (U) - (5 cos (U 



2sin2 ( 1(5 



(75) 

(76) 
(77) 



We will derive equality (75) from equality (77). It follows from (73) 
that 

1 - cos^ (U) 

C2 = (78) 

1 + ^2-20; cos I IS 
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Solving the corresponding quadratic equation for cos yj^j ^rid tak- 
ing into account that cos (\^^ < 0/ we get 



cos 



] =ciC^- x/cc^C* - C2 - a2C2 + 1. 



Therefore, 



cos ( ^(5) =aC2--^l-C2^1-a2c2_ (79) 



Next, using (78) and (79), we see that equality (75) is equivalent to 
the following: 



Xl - Xq 
Vi 



l + fl;2_2acos (^^S 



\ sin2 (1 



+ 



l + a2-2acosfM) S{1 + oc^ - la cos (is)) 
. . ^ ^ - 1 + , , ^ ^ ■ (80) 

2^ 



\ sin2(^icf) 2sin2(^i^ 



Making simplifications in (80), we can show that equality (80) re- 
duces to equality (77). 
This completes the proof of Lemma 13. 

Lemma 14. The S-formula for the Heston distance in the case, where the 
■points {xq,vq) G % and {x\,v\) G % are S-far, follows from part (ii) of 
Theorem 3. 

Proof. In the proof of Lemma 14, we will use the notation intro- 
duced in the proof of Lemma 13. Suppose the points {xq, vq) G H 
and {xi,vi) G Ti are ^-far, and assume that part (ii) of Theorem 3 
is valid. By Lemma 12, {xo,vq) and {xi,Vi) are C-far, and it follows 
from part (ii) of Theorem 3 that 

arccos(C*^/uo) + arccos(C*v^i)) 
dH{ixo,Vo),{xi,vi)) =2 ^ ^-^^ ^ i^. 

Taking into account (71), (73), and (79), we see that 

-arccos(aC) + arccos(C) 



%((^0/^o)/ {xi,vi)) = 2v^ 

= 



C* 

arccos(aC2 - Vl - C2 Vl - a2c2) _ 



C* c 
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Next, using (73), we obtain 



S^JviJ 1 + (X^ — 2(X cos ^2^^ 



sin (^^S 



Finally, recalling that a. = we see that formula (32) holds. 
This completes the proof of Lemma 14. 



5. The limiting cumulant generating function for the 

Grushin model 

It is not hard to see using the equation in (28) that the Laplace 
transform pp((xo/i/o)/ (^/i/l)) iri the variable Xi of the Grushin tran- 
sition density pj satisfies the following heat equation with quadratic 
potential: 

with the initial condition given by ((xo/J/o)/ = ^xo{^i)^^'^^°- 

The fundamental solution for such a heat equation is well-known 
(see, e.g., [4], Theorem 10.3). Using this fundamental solution, we 
get 



w 



71 sm 



(f) 



, ^ [-X0 sin (f ) + {y\ + yl) cos (f ) - lyry^] 

sin(f) 



The function pf given by (81) has a removable singularity at w = 0. 
The analyticity strip for pf is given by — ^ < ?R(w) < It is not 
clear yet what happens outside this strip. 
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We will next compute the Laplace transform of the function pf in 
the variable i/i, using formula (81). We have 



00 ("OO 



00 J —00 



e'''^+^y^pf{{xo,yo),{xiAJi))dyidxr 

z(;(-xosin(f)+y2cos(f)) 



w 



71 sm 



e^^i exp 



(?) 



exp 



sm 



(f) 



IV 



[y^os(f)-2yiyo] | 

— ^i;^?^ r^^ 



Let us next replace w by j and /3 by j in the previous equality. This 
gives the following: 

/CO l-CO ( ^ 

/ exp <^ -{5xx + 7yi) \ pp((xo,yo), (xi,yi))dyidxi 
-oo J — 00 I r I 



^ nt sin ^1 



exp 



^ (^-xq sin j + y^ cos 



f sin U 




exp 



exp 





COS 


(i) 


- 


2yiyo 


f sin ( 







dyi. (82) 



The new restrictions on the parameters are —2n<S<2n and 7 G 
R. 

Denote the expression on the left-hand side of (82) by }t{xo, S, yo/ 7) • 
Assume —7T<S<n and 7 G R. After lengthy but straightforward 
computations, we obtain 

1 

/i(^0/<5,yo,7) = 



cos 



X exp ^ y 



I {AS^yl + 72) sin (|) + AxqS^ cos (|) + 4yo7^ 



43 cos 



For every pair (S, 7) G IR^ ^^d t > 0, put 

At (^,7) = ilog/t(xo,^,yo/7) 

and 

A{3,j) = limAi(^,7) 



(83) 

(84) 
(85) 
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if the limit in (85) exists. Note that the functions and A depend on 
Xq and 1/0- The function A is called the limiting cumulant generating 
function associated with the Grushin model. 

Let us assume that — zr < ^ < zr and 7 G R. Then, taking the 
logarithm of the integral on the left-hand side of (83), multiplying 
by t, and sending t to infinity, we obtain 



(4^2y2 ^ ^2^ U\ + 4xoS^ COS (1) + 4yojS 

= — 7-x — 

4^ cos (I) 



(86) 



Let us next assume that ^ G (—271, — tt] IJ [n,2n) and 7 G R. Then, 
using (82), we see that Af((5, 7) = 00 for every t > 0. Here we take 
into account that under the restrictions imposed above, 

^cos (^0 sin-1 (^0 <0. 

It follows from (85) that A(^, 7) = 00 for all S G {-2n, -n]\J [tt, 2n) 
and 7 G R. Now, using Holder's inequality, we see that A(^, 7) = 00 
for all S G (—00, — tt] U [tt, 00). 

The limiting cumulant generating function A is defined every- 
where and convex on R^. This function is finite on the set D = 
(— TT, tt) X R and identically infinite outside this set. Moreover, the 
function A is continuous on the set [R\ { — tt} \ {tt}] x R. It is also 
continuous on the lines 3 = n and 3 = —n with the exception of the 
points Pi = (tt, — 2711/0 ) and P2 = (— tt, —2nyo). More precisely, we 
have A (Pi) = A (P2) = 00. In addition, 

liminf A(^,7) = tixq (87) 



and 



liminf A(3,y) = -nxn. (88) 



Hence, the limiting cumulant generating function A is lower semi- 
continuous everywhere in R^, except at the points Pi and P2. 

Using the definition of the function A, we see that for all (^, 7) G 



4i/^^2 _^ sin ^) + ^'^ (5 - sin ^) + 4^2^yQ 



+ X0 (89) 
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and 



2^ COS ( f 



(90) 



Hence, the function A is continuously differentiable on the set D. 
However, this function is not steep (the definition of the steepness of 
a function is given in [10]), Definition 2.3.5). We will next prove the 
previous statement. 

It follows from (90) that 1 1 V A (c^, 7) 1 1 ^00 provided that {S, 7) G 
D and {S, 7) — > {Sq, 70) with either Sq = n and 70 7^ —2nyo, or Sq = 
— 71 and 7o 7^ —2nyo. The behavior of the gradient at the exceptional 
points Pi = (tt, —Inyo) ^rid P2 = ("^r, —Inyo) can be described us- 

aA 



ing (89), (90), and I'Hopital's rule. We have 



and 



aA 
a7 



.2yo 

71 



dS 



27T 



71 



as (c^,7) G D and (c^, 7) ^ P\. Therefore, 



4+^0 



|VA ((5,7) 



\ 




4) + Xo 



as (^,7) G D and ((5, 7) — > Pi. Similarly, 



aA 

d5 



^(4 

271 



71 



dA 2yo 

a7 71 



and 



|VA ((5,7) 



\ 




71- 



■) + Xo 



as ((5, 7) G D and {5, 7) — ^ P2. Therefore the steepness condition 
for the function A is satisfied everywhere on the boundary of the set 
D(A), with the exception of the points Pi and P2. 

Remark 15. The absence of the lower semi-continuity and of the 
steepness property for the function A does not allow us to use the 
Gartner-Ellis theorem (see Theorem 2.3.6 in [10]) to establish the large 
deviation principle for the Grushin model. It will be shown below 
that this principle is valid only in a special regime (see Theorem 22). 
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6. Critical points and the Legendre-Fenchel transform 

Let us consider the Legendre-Fenchel transform A* of the limiting 
cumulant generating function A. It is given by 

A* (xo, 1/0/ ^1/1/1 ) = sup [{xi,5) + (i/i,7> - A(^,7)]. 
It is clear that 

A*(xo,i/0/^i/!/i) = sup + (1/1,7) - A((5,7)]. (91) 

-7r<(5<7T,76R 

A critical point (c^, 7) = ((J*,7*), —n<5* < n, for the function 
(xi, 5) + (1/1, 7) — A(^, 7) satisfies the following system of equations: 

aA 

It follows from (92) and (90) that 

25 (yi cos (f ) - yo) 
7 = — ^ -H: -■ (93) 



sm I 2 



Moreover, (92) and (89) imply that S* is a solution to the following 
equation: 

4y2^2 ^ sin(^)) + 7^ - sin(^)) + AS^jyo sin 



8(52 cos2 (^1^ 

= xi - xo. (94) 

The second component 7* of the critical point ((5*, 7*) can be found 
by plugging the solution S* to the equation (94) into (93). This gives 

2S* (yicos(^)-yo) 
7* = ^ ^ -■ (95) 



sm 



Our next goal is to simplify the equation in (94) by taking into ac- 
count (94) and (95). Tedious, but rather straightforward calculations 
show that (94) can be rewritten in the following form: 

f{S) = f{5, yo, yi) = xi - xo, (96) 

where 

/(^) = \Mm5) (97) 
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with 



A{S) = A{S,yo,yi) = sin 



-2 



( 



2 



S 



) 



(98) 



and 




(99) 



It is assumed in (96) that —7t<S<7t. The value of the function / at 

1 . _ 

(5 = is computed as follows: /(O) = - \im A{S)B{S) = 0. 

Remark 16. It follows from (97) - (99) that the function / is defined 
on the interval (—271, 2 tt). This fact will be used below. 

Remark 17. The properties of the functions /, A, and B are simi- 
lar to those of the functions /, A, and B defined by (52), (53), and 
(54), respectively. The previous statement follows from the equalities 

/(^/yo,yi) = fi^.ylvl), M^'VO'Vi) = M^'VlyD' ■B(^,yo,yi) = 

H^'Vo'Vi)' and S*{{xo,yo), = 5{{xo,yo), (xi,i/i)). 

Let us fix xq and i/o- h follows from Lemma 2.3.9 in [10] that the 
Legendre-Fenchel transform A* of A is a good rate function. Explicit 
formulas for the function A were found in Section 5. In the present 
section, we compute the function A*. A simple analysis of formula 
(91) defining the function A* shows that to compute the supremum 
in (91) one has to take into account the input of the critical point 
(^*,7*), the boundary of the strip where the moment generating 
function is finite, and the boundary at infinity. Since A{5,'y) — > oo 
as 7 ^ 00 or 7 —> —00, the input of the boundary at infinity can be 
ignored. Using formulas (87) and (88), we see that the the input of 
the exceptional points Pi and P2, more precisely, of sequences con- 
verging to those points, is given by the following expression: 



Note that the number in (100) is positive if and only if |xo — xi| > 

2yiyo- 

Next, suppose S* ^ (— tt, zr). This means that 



^1 = 7r(|xo - ^il -2yiyo)- 



(100) 




(101) 
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In this case, there is no critical point inside the fundamental strip, 
and hence A* is given by 



A*(xo,i/0/^i/J/i) = 7r(|xi - xol -2i/ii/o) = ^i- 



(102) 



On the other hand, if ^* G (— tt, n) (this happens if the opposite 
inequality to the inequality in (101) holds), then the input that the 
critical point (^*,7*) brings to the computation of the supremum in 
the formula for A* is given by the following expression: 



(4(^*)2y2 + (7*)2) sin(^)+47*^*yo 



43* cos 



(103) 



Replacing xi — xq in formula (103) by the expression on the left-hand 
side of formula (96) and making simplifications, we obtain 



*\2 



Vi + yl- 2yiyo cos 



(104) 



Therefore the condition | xi — xq I < 2yiyo implies the equality 

A*(xo,yo,^i,yi) = ^2, (105) 
and the condition 2yiyo < |xi — xo| < 5 [tt (y^ + yg) + 4yiyo] gives 



A*(xo,yo,Xi,yi) = max{Ri,R2} 



(106) 



Our next goal is to compare all the inputs discussed above. We 
will next show that R2 always dominates Ri. 

Lemma 18. For all (xo,yo) £ IR^ i^nd (xi,yi) G R^, the following in- 
equality holds: R\ < R2- In addition, if 



^[yl + yo) +4yiyo 



< Fi - xo\, 



(107) 



then Ri < R2. 



Proof Taking into accound the definition of Ri and R2 (see (100) 
and (103), respectively), replacing the expression xi — xq in formula 
(103) by the expression on the left-hand side of formula (96), and 
simplifying, we see that the inequality Ri < R2 can de derived from 
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the inequality 



ln5* cos 



+ yiyo 
> (vl + vt) + yiyo 



+ 4/1 sin 



*\2 



COS 



, + 471 sin , 
2 J \ 2 



(108) 



It is easy to see that with no loss of generality we may assume that 
5* > 0. We will prove that 

(^*)2 + 7rsin((5*) > nS* 



(109) 



and 



nS* cos 



+ 27rsin^ 



> (^*)2cos 



+ 271 sin 



5* 
~2 



(110) 



for all < ^* < 71. Moreover, it will be shown that the condition 
7r < ^* < 271 implies strict inequalities in (109) and (110). Note that 
condition (107) is equivalent to the condition n < 5* < 2n. It is clear 
that Lemma 18 follows from the inequalities formulated above. 

Suppose first that < ^* < 7r. We will next establish (109). The 
fact that the inequality in (110) is equivalent to the inequality 



d*{n - cos ( ^ 1 > 27rsin ( 



sm I — 



(111) 



will be used in the proof. 

Let us assume that (109) holds for all < ^* < f , and let f < 
(5* < n. Then d = n — d"" satisfies (109), and it is easy to see that ^* 
also satisfies (109). It follows that it suffices to assume < cJ* < ^. 

Using the Taylor series, we see that for < cJ* < ^, 



sin(^* 



> 1 



cJ* - 6 
Hence (109) can be derived from the inequality 



n 1 



> n 



The previous inequality is equivalent to 71(5* < 6, which is of course 
correct. This establishes (109). 
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Our next goal is to prove (111). Let us first assume < S* < j 
Using the Taylor series, we obtain 



cos ( y 1 > 1 



g and sin I — I > — 



48 



(112) 



Therefore, (111) can be derived from the following inequality: 



*\2 



8 



> 1 



*\3 



2 ^ 48 



The previous inequality is equivalent to 

6{S*)^ + 2471 > 48 + 671^ + n{S*f. 
It is not hard to see that (113) can be rewritten as follows: 



(6-71) 



3n 

6—71 



+ 



2471 - 48 



9n^ 



6—71 



Since S* < f , (114) follows from 



(6-71) 



371 71 

6-71 ~ 2 



+ 



2471 - 48 

6—71 



(6-71) 



971^ 

(6-7r) 



> 0. 



(113) 



(114) 



> 0. 



(115) 



Now, it is clear that (111) holds for all < S* < ^, since (115) can be 
easily checked. 

We will next show that (111) also holds under the condition ^ < 
3* < 71. The following estimate will be needed in the sequel. For all 
J with < <5 < f. 



d{n — 5) sin ( - I > 27rcos 



1 — cos 



(116) 



Dividing the both sides of (116) by cos ^1 j and using the inequality 

tanx > X for < x < ^, we see that (116) can be derived from the 
inequality 

(117) 



5\7i-5) > 4:71 (^1 - cos 

Next, using the first inequality in (112), we see that (117) follows 
from the inequality 

{n-S)>^^ (118) 

which clearly holds because < S < j. 

Suppose that ^ < S* < 7Z and denote S = 7Z-S*. Then < J < f , 
and hence inequality (118) holds for S. It is not hard to see that (118) 
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for S is equivalent to (111) for S*. This completes the proof of estimate 
(110). 

It follows from (109) and (110) that esimate (108) holds. It has 
already been mentioned that (108) implies the inequality Ri < R2. 
Therefore, part of Lemma 18 in the case where — zr < ^* < zr is 
valid. 

Now let TT < (5* < 2tz. We will first establish that the strict inequal- 
ity in (109) holds. It is clear that the function pi{u) on the left-hand 
side of (109) and the function p2{u) on the right-hand side equal tt^ 
at u = TT. Moreover, p[{u) = 2u + ncosu and P2{u) = tt. It is easy 
to see that p\{u) > P2{u) for all n < u < In. It follows that the 
strict inequality in (109) holds when n < S* < In. The proof of the 
strict inequality in (110) under the same restriction is similar. Here 
the functions p\ and p^ equal 2/1 at w = n. Moreover 

, . , u n u u u 

Piiu) = 71 cos usin — |-27rsm — cos — 

riv / 2 2 2 2 2 



and 



/ / \ U 1 J . u ^ u 

Pj[u) = TTcos — — -u sm — + 2m cos — . 
ri\ J 2 2 2 2 



It is not hard to see that p'i(u) > P2{u) for all n < u < In. This 
implies the strict inequality in (110) in the case where n < 5* < In. 
The proof of Lemma 18 is thus completed. 



Theorem 19, Under the condition 

1 



Fi - ^ol < 



2 L 



7^ (yi+yo) +4yiyo 



(119) 



the Legendre-Fenchel transform A* of the limiting cumulant generating 
function A in the Grushin model is given by the following formula: 



A*{xo,yo, xi,yi) 



*\2 



2 sin 



lis 



yl + yl-^ym COS ( — 



(120) 



On the other hand, if 



1 
2 L 



then we have 



(121) 



A*{xQ,yQ,xi,yi) = 7r(|xi - xqI - 2yiyo)- 



(122) 
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Remark 20. Recall that 3* in formula (120) is the unique solution to 
the equation 



1 

- CSC - 



2 

Xi - Xq 



(y? + yo) - sin(^)) - 2yiyo cos (^0 - 2 sin (^0 ^ 



The solution to this equation satisfies — zr < (5* < tt if condition (119) 
holds. 

Proof of Theorem 19. Formula (122) in Theorem 19 has already been 
established (see (102)). Formula (120) can be derived from (104), 
(105), (106), and Lemma 18. 

Remark 21. Let us consider a special case where Xi = xq. Then 
3* = and therefore (120) gives 

A*{xo,yo,xi,yi) = 2{yi - y^f. (123) 



7. A PARTIAL LARGE DEVIATION PRINCIPLE FOR THE GRUSHIN 

MODEL 

Le us recall that in the present paper we denoted by A and A* 
the limiting cumulant generating function in the Grushin model and 
the Legendre-Fenchel transform of A, respectively. We will prove 
below that for any initial point A = (xo/yo) ^ IR^ the large deviation 
principle holds for the Grushin model in a certain open subset M of 
the plane R^. The set M consists of all the points which are ^-close 
to A. In the formulation of the next theorem, the symbols B° and B 
stand for the interior and the closure of the set B, respectively. 

Theorem 22. Let xq ^ IR md yo G ]R he fixed, and consider the open set 
in defined by 



M 



xi,yi) GR^: |xi-xo| < tt yf + y^ + 4yiyo ^.(124) 



Then the large deviation principle with the rate function A* holds on the 
set M. More precisely, for any Borel subset B ofM, 

- inf A*(xo,yo/^i/yi) < liminf flogP/^(xo/yo/B) 

(xi,yi)€B° t^O L 



< lim sup t log Pf ( Xq, yo/ B ) 



<- inf _A*(xo,yo,^i,yi)- (125) 

(xi,yi)6B 
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In addition, if the set B is such that B° = B, then 



lim 



tlogPt{xo,yo,B) = - inf A*{xo,yo,Xi,yi). (126) 

J (xi,yi)€B 



Proof. We have shown in Section 6 that for all {xi,yi) G M, the 
following equality holds: VA(^*,7*) = {xi,yi). Therefore, Lemma 
2.3.9 in [10] implies that any point (xi,i/i) G M is an exposed point 
of A* with ((5*, 7*) being the exposing hyperplane for {xi,yi) (see 
[10] for the definition of exposed points and exposing hyperplanes). 
In other words, the set M consists entirely of exposed points. Now it 
is not hard to see that all the conditions in the Gartner-Ellis theorem 
(see Theorem 2.3.6 in [10]) hold. Applying this theorem and taking 
into account the continuity of the function A*, we establish Theorem 
22. 

8. The distance formula 

Our goal in this section is to complete the proof of Theorem 3. 
Since we have already proved this theorem in the far ^-regime, it 
remains to establish it for the pairs of points {xq,Vq) G 1-1° and 
{xi,Vi) G H° in the close ^-regime. 

We will prove the following assetrion: 

Theorem 23. If condition (119) holds, then 



A*{xo,yo,Xi,yi), (127) 



> A*{xo,yo,xi,yi). (128) 



4((^o^yo)^(^i^yi)) 

2 

while if condition (121) is valid, then 

2 

Remark 24. It is clear that Theorem 3 in the close ^-regime follows 
from (30), (120) and (127). Inequality (128) is interesting because it is 
often expected that for the distances associated with stochastic mod- 
els, the function jd^ coincides with the Legendre-Fenchel transform 
A* of the limiting cumulant generating function A. Inequality (128) 
shows that this is not the case for the Grushin model when the points 
are in the far ^-regime. 

Proof of Theorem 23. Since the Grushin operator (27) is hypoelliptic, 
it follows from the results obtained by Leandre in [ j], [Z4] (see also 
[25], [26]) that Varadhan's equality, that is, the equality 

lim(nogpP((.„,y„),(xi,yi))l = - 4((^o,yo),(xi,y.)) (1^9) 
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holds uniformly on compact subsets of (Varadhan's results can 
be found in [29, 30]). 

Fix a point (xo,i/o) £ IR^ and consider all the points 
such that condition (119) holds. Let B = Be(xi,i/i), where B£(xi,i/i) 
is the disk of radius e in centered at (xi,i/i) and such that 

Beixi,yi) C M 

(recall that M is defined in (124)). It follows from (126) that 



lim [t log Pt (xo, yo/ Be(^i/!/i))] 
= ~. ^ A*(xo,yo,:^,y) 

{x,lj)eBe{xi,lJi) 



(130) 



for all small enough e > 0. Using the mean value theorem for inte- 
grals, we can prove that 



tlogP^ {xo,yo, Be{xi,yi)) + 



< 



ilog(7re^) 



+ sup 

{x,y)eBe{xi,yi) 

+ sup 

{x,y)eBe{xi,yi) 



4((^o,yo),(^,y)) 



t log ft {xQ,xi,x,y) + 



4((^o,yo), {xi,yi)) 4 ((^o,yo), {x,y)) 



(131) 



It follows from (131) and (129) that 



lim sup 



,1 dG/ n t ^^ , 4((^o,yo),(:^i,yi)) 
ilogPf^ {xo,yo,Be{xi,yi)) + 



< sup 

(x,y)€Be(xi,yi) 



4((^o/yo)/ (^1/yi)) 4 ((^o,yo), (^,y)) 



(132) 



Taking the limit as £ in (132) and using the continuity of the 
function cfi (the continuity follows from (31)), we get 



lim lim sup 



vGf u ( ^^ , 4((^0/yo)/(^i/yi)) 
ilogPf^ (xo,yo, Be{xi,yi)) + 



0. 



(133) 
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It is not hard to see that 

4((^o/yo)/(^i/yi)) 



< 



2 



tlogP^ {xQ,yo, Be{xi,yi)) + 



A*(xo,!/0/ ^i/J/i) 

4((^o,yo),(^i,yi)) 



+ 



+ 



tlogPt {xo,yo,Be{xi,yi)) + inf A* (xo,yo/ ^/y) 

(x,y)6B£(xi,yi) 



inf A* {xo,yo,x,y) - A*{xo,yo,xi,yi] 

(x,ij)eBe(xi,iji) 



(134) 



Taking the limit as e ^> in (134) and using (130), (133), and the 
continuity of the function A*, we obtain formula (127) for all pairs 
of points (xo,yo), {xi,yi) satisfying condition (119). In addition, it 
is not hard to see that formula (128) for all pairs of points (xo,yo)/ 
satisfying condition (121) follows from (122), (100), (104), the 
second part of Lemma 18, and from Theorem 3 in the far ^-regime. 
Note that we have already established Theorem 3 for points which 
are S-iar. 

This completes the proof of Theorem 23. 
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